arXiv:1505.02139vl [astro-ph.SR] 8 May 2015 


Mon. Not. R. Astron. Soc. 000, [T]-?? (2000) Printed 11 May 2015 (MN style file v2.2) 


First detections of 610 MHz radio emission from hot magnetic stars 

P. Chandra^*, G. A. Wade^, J. O. Sundqvist^, D. Oberoi^ J. H. Grunhut"^, A. ud-Doula 
V. Petit^, D. H. Cohen^, M. E. Oksala^, and A. David-Uraz^’^ 

^National Centre for Radio Astrophysics, Tata Institute of Fundamental Research, P.O. Box 3, Pune 411007, India 
^Department of Physics, Royal Military College of Canada, PO Box 17000, Station Forces, Kingston, Ontario K7K 7B4, Canada 
^Department of Physics &• Astronomy, University of Delaware, Newark, DE19716, USA 
‘^European Southern Observatory, Karl-Schwarzschild-Str. 2, D-85748 Garching, Germany 
^Penn State Worthington Scranton, 120 Ridge View Drive, Dunmore, PA 18512, USA 

^Dept. of Physics &- Space Sciences, Florida Institute of Technology, Olin Physical Science, 346, Melbourne, FL 32901, USA 
Department of Physics and Astronomy, Swarthmore College, Swarthmore, PA 19081, USA 

^LESIA, Observatoire de Paris, CNRS UMR 8109, UPMC, Universite Paris Diderot, 5 place Jules Janssen, 92190 Meudon, France 
'^Department of Physics, Engineering Physics and Astronomy, Queen’s University, 99 University Avenue, Kingston, ON K7L 3N6, Canada 


Submitted 8th May 2015 


ABSTRACT 

We have carried out a study of radio emission from a small sample of magnetic O- and B- 
type stars using the Giant Metrewave Radio Telescope, with the goal of investigating their 
magnetospheres at low frequencies. These are the lowest frequency radio measurements ever 
obtained of hot magnetic stars. The observations were taken at random rotational phases in 
the 1390 and the 610 MHz bands. Out of the 8 stars, we detect five B-type stars in both 
the 1390 and the 610 MHz bands. The O-type stars were observed only in the 1390 MHz 
band, and no detections were obtained. We explain this result as a consequence of free-free 
absorption by the free-flowing stellar wind exterior to the closed magnetosphere. We also 
study the variability of individual stars. One star - HD 133880 - exhibits remarkably strong and 
rapid variability of its low frequency flux density. We discuss the possi bility of this emission 
being coherent emission as reported for CU Vir bv lTrigilio et al.l (l2000l) . 

Key words: radiation mechanismsmon-thermal — stars; massive — stars: individual: HD 
133880, etc. — stars: magnetic held — radio continuum: stars 


1 INTRODUCTION 

The last decade has witnessed the identification and elaboration of 
an important population of hot O - and B-type stars hosting strong, 
organ i zed magnetic fields (e.g. DonatieGD 1 20021 : lOksala et al.l 


bol d: [Grunhut et^l2012l : IWade et alj 201 j) . Observational and 


theoretical studies of the interaction of their intense radiation- 


driven winds with their magnetic fields (e.g. 

ud-Doula & Owockil 

I 2 OO 2 I; iTownsend et alj|2005l; ISundavist et alj 

2 OI 2 I) have revealed 


that the stellar wind properties of magnetic hot stars are modified 
in important ways as compared to the winds of non-magnetic stars, 
introducing significant changes across the entire electromagnetic 
spectrum. 

The presence of an organized magnetic field at the surface of 
a hot star leads to channelling and confinement of its outflowing 
stellar wind, creating a magnetosphere, which can radiate in various 
wavebands dAndre et al.lll988l : lLinskv 

In regions close to the star, the magnetic pressure dominates 
and causes the wind to follow dipolar field lines. At greater dis¬ 
tances from the stellar surface, the wind kinetic energy density 
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exceeds the magnetic pressure due to the stronger radial decline 
of magnetic field energy density than of wind kinetic energy den¬ 
sity. The radius at which the wind kinetic energy density becomes 
equal to the magnetic press ure is defined as the Alfven radius (e.g. 
lud-Doula & Owockill2002l) . The region below the Alfven radius, 
i.e. interior to the lar gest closed magnet ic loop, is defined as the 
inner magnetosphere ITrigilio et al.l2004l). which is the site of gen¬ 
eration of X-ray and Ha emission (e.g. Babel & MontmerI3l 19971 : 


iGagne et al.l2005l : iHowarth et alj20(^ : ISundgvist et alj201^ . Be- 

yond the Alfven radius, the wind opens the magnetic field lines, 
generating a current sheet in the magnetic equatorial plane. This 
region is defined as the middle magnetosphere. At greater distances 
from the star, the outer magnetospheric region is characterized by 
a radial magnetic topology following wind streamlines. 

Radio emission from non-magnetic hot stars is expected to be 
thermal free-free emission from the ionized wind in the circumstel - 
lar environment dPanagia & Fellilflg^ : IWright & Barlo^ll975h . 
However, in the presence of magnetic fields electrons can be accel¬ 
erated to relativistic energies, either by reconnection near the c ur- 
rent sheet in the middle magnetosphere dUsov & Melroseiri992l) or 

[ enerated in strong, large-scale shocks in the inner mag netosphere 
Babel & Montmerljl 19971 : lud-Doula et ahlbOOd 12008 ). These en- 
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ergetic electrons can give rise to gyrosynchrotron radio emission 
JOrake et al.ll 19871 : iLinskv et~^l 19921 : iTrigilio et alj2004) . 


Various surveys of cooler magnetic A- and B-type stars have 
been carried ou t at 5-8 GHz radio frequencies, result i ng in ^ 25% 
detec ti ons (e.g. Abbott et^ll985l : lOrake et ^1 19871 : iLinskv et~n] 
Il992h . ILinskv et alJ (Il992h interpreted the radio emission in these 
stars as gyrosynchrotron in nature, arising from regions located at 
10-20 Rt, with higher frequency emission originating from plasma 
located closer to the star, and the lower frequency emission from 
material further away. 

iLeone & Umanal j 19931) monitored the 5 GHz radio emission 
of two magnetic B-type stars, HD 37017 and a Ori E, and ob¬ 
served variability in accordance with their respective rotational pe¬ 
riods. The coincidence of radio maxima with the extrema of the 
longitudinal magnetic field in these stars led those authors to sug¬ 
gest that radio emitting regions are located above the magnetic 
poles. Radio variability was also measured from the magnetic B8p 
star HD 133880 (HR Lup). The flux density and circular polariza- 
tion were observed to vary accordin g to the stellar rotation period 
dLim et alJI 19961 : iBailev et^ 20121). and was also interpreted as 
gyrosynchrotron emission l lBailev 'rtalj2012h . 

iLeone et"^ ( |2004|) argued that as per the standard radio emis¬ 
sion gyrosynchrotron models, emission at a given frequency should 
be emitted in a well-localised torus above the magnetic poles. They 
examined the radio spectral energy distributions (SEDs) of five 
magnetic B-type stars in the range from 1 GHz to 87.7-GHz. They 
concluded that the observed slopes of the radio spectra and the ab¬ 
sence of millimetre emission are not generally compatible with this 
model, and suggested a cut-off frequency of the radio emission. 


ITrigilio et al.l ( |2000|) detected rapid, intense, narrow-band and 
highly polarized radio bursts from the late B-type star CU Vir, 
which they hypothesised to be Electron Cyclotron Maser Emission 
(ECME) at 1.4 GHz. Such emission is expected to occur principally 
at low frequencies (< 1.5 GHz, ITrigilio et al.l I 2 OO 0 I : iLeto et al.l 
I2OI2I) . However, to date, few such low frequency observations have 
been obtained, and gene rally they correspond to null results (e.g. 
ICeorge & Stevensll2()12l) . 


Clearly, a variety of diverse investigations of the radio emis¬ 
sion of magnetic A- and B-type stars have been carried out, lead¬ 
ing to an array of disconnected results. However, no similar studies 
have yet taken place for hotter magnetic O- and B-type stars. In an 
attempt to extend and homogenize the study of the physics of radio 
magnetospheres of hot stars, we have initiated a systematic survey 
of the radio emission properties of the known magnetic B- and O- 
type stars. Our survey includes observations at both high frequency 
(currently being obtained with the VLA) and at low frequency (cur¬ 
rently being obtained with the GMRT). 

In this first paper we describe a pilot study of 8 established 
magnetic O- and B-type stars, observed with the GMRT in the 
1390 MHz (20 cm) and 610 MHz (50 cm) bands, to characterize 
the radio emission properties of their magnetospheres at very low 
frequencies, study their radio SEDs to understand magnetospheric 
properties and probe mass loss from the star in the case of dominant 
free-free absorption. 

In ^we describe the observations and data analysis. In ^ 
we carry out a detailed analysis and look for variability of the de¬ 
tected stars within the observing period, and in g3.6l we discuss the 
stars for which no radio emission was detected. In 0we discuss 
our results and the potential of low frequency radio observations 
for providing new constraints on the physics of hot star magneto¬ 
spheres. 


2 OBSERVATIONS 

The GMRT observations of 8 magnetic O- and B-type stars were 
taken between 2013 October 24 to 2014 January 24 under cycle 
25 in the 1390 (L-band) and 610 MHz bands. Table[T]gives details 
of the observed stars. Out of these, 5 B-type stars were observed 
in both the 1390 and 610 MHz bands, while the 3 O-type stars 
were observed only in the 1390 MHz band. Each observation was 
of approximately 3 hours in duration and a total of 26-29 good an¬ 
tennae were used. Eor all observations, visibilities were recorded 
for two circular polarizations (RR and LL) with a bandwidth 32 
MHz, divided into 256 frequency channels, and 16 s integration 
time. Calibrator sources were used to remove the effect of varia¬ 
tion of the instrumental and other non-astronomical factors in the 
measurements. 3C48, 3C286 and 3C147 were used as flux cali¬ 
brators in various observations. Elux calibrators were observed for 
10-15 minutes, either in the beginning or towards the end of each 
observation. Elux calibrators were also used for bandpass calibra¬ 
tion. Phase calibrators were chosen from the VLA calibrator man¬ 
ual such that they are located within 15 degrees of the target star. 
They were observed more frequently. This is important not only for 
the tracking of instrumental phase and gain drifts and atmospheric 
and ionospheric gain and phase variations but also for monitoring 
the quality and sensitivity of the data and for spotting occasional 
gain and phase jumps. In the 1390 MHz band, the phase calibrator 
observations were made approximately every 25-30 minutes for 
5 minutes duration, whereas in the 610 MHz band, the 6-minute 
phase calibrator scans were obtained every 35^0 minutes. 

A ELAGging and CALibration (ELAGCAL) software pipeline 
devel oped for automatic flaggi ng and calibration of the GMRT 
data dPrasad & Chengaiu3l2012l) was used to flag and calibrate the 
GMRT data. In 70% of the cases, the results were satisfactory and 
only minor additional flagging was needed after the ELAGCAL. 
However, in some cases, the data needed more significant flagging 
and recalibration. In such cases the flagging was done manually and 
a flag file was created. The ELAGCAL was then run again using 
the manual flag file. The calibrated data were inspected to deter¬ 
mine the presence of any spurious signature. The above step was 
repeated until satisfactory calibration was achieved. The flagged 
and calibrated visibility data were used to make continuum images 
using the standard tasks in the Astronomical Image Processing Sys¬ 
tem (AIPS). To avoid bandwidth smearing in the continuum im¬ 
age, the total bandwidth was divided into 6 sub-bands in the 1390 
MHz observations, and in 10 sub-bands in the 610 MHz observa¬ 
tions. Because of the large field of view (EoV) of the GMRT, the 
three-dimensional (3D) imaging feature in the AIPS task TMAGR’ 
was used in which the entire FoV is divided into multiple subfields 
(facets) and each of which was imaged separately. For the 1390 
MHz image, the total FoV was divided in 19 subfields, and in the 
610 MHz band the FoV was divided into 37 subfields. The pres¬ 
ence of a large number of bright sources in the FoV of various 
stars allowed us to carry out self-calibration to improve the com¬ 
plex gains. This reduces the errors from temporal variations in the 
system gain, and spatial and temporal variations in the ionospheric 
properties. After 3 rounds of phase self-calibration, the clean com¬ 
ponent model was subtracted from the UV data to identify the 
residual bad data. Some more flagging was performed and 2 more 
rounds of phase self-calibration were run. A final round of ampli¬ 
tude and phase self-calibration was also performed in each dataset. 
Since the phase variations occur on time-scales of a few minutes, 
the time interval for the phase self-calibration was chosen to be 1 
minute, and 5 minutes for the amplitude and phase self-calibration. 
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Table 1. Details (spectral type, eflFective temperature, luminosity, radius, mass, rotational period, polar strength of the magnetic dipole at the stellar surface, 
Kepler co-rotation radius, Alfven radius, mass loss rate and wind terminal velocity) of magnetic O and B stars observed with the GMRT. Rf^ is the radius of 
the free-free absorption photosphere at 1390 MHz, discussed in Sect. m 


Star 

Spectral 

Teff 

logT* 

R* 

M* 

Period 

fip 

Rk 

Ra 

log(M) 

Vco 

Rff 

Rff/RA 


Type 

(kK) 

(To) 

(Re) 

(Me) 

(d) 

(kG) 

(R.) 

(R.) 

(Me yr-’) 

(km s“’) 

(R.) 


HD215441 

B8-9p 

15 ±2 

1.8 

2.6 

3.5 

9.49 

34 

- 

- 

- 

- 

- 

- 

HD37479 

B2 Vp 

23 ±2’ 

3.6’ 

3.9’ 

5.0’ 

1.19’ 

9.6” 

2.1” 

31” 

-9.774° 

1794° 

6.5 

0.2 

HD37017 

B2 Vp 

21 ±2'’ 

3.4'’ 

3.9'’ 

7.2'’ 

0.90'’ 

>6.0” 

1.9° 

> 18° 

-9.072° 

1102° 

27.6 

< 1.5 

HD36485 

B3 Vp 

20 + 2’ 

3.5’ 

4.5’ 

7.1’ 

1.48’ 

10’ 

2.4° 

24° 

-8.969° 

1012° 

30.4 

1.3 

HD133880 

B8IVp 

13 + 1” 

2.0” 

2.0” 

3 2” 

0.88” 

19” 

2.2” 

60” 

-ll.O” 

750” 

4.6 

0.1 

HD37022 

07 Vp 

39 ± 1^ 

5.3« 

9.9'’ 

45'’ 

15.42’° 

1.1” 

9.4° 

2.4° 

-6.399® 

3225° 

236.3 

98.5 

HD57682 

09 V 

34± 1’2 

4.8” 

7.0” 

17” 

63.57’” 

1.7’” 

24° 

3.7° 

-7.079® 

2395° 

153.7 

41.5 

NGC 1624-2 

06.5f?cp 

35 ±2’'’ 

5.1’'’ 

97’4 

34’'’ 

158.0’'’ 

> 20’'’ 

41’4 

> 11’'’ 

-6.786° 

2890° 

151.2 

< 13.7 


’Hunger et al. 1989, ^Townsend et al. 2010, ^Oksala et al. 2012, '’Bolton et al. 1998, ^Bohlender et al. 1987, ®Petit et al. 2013, ^Leone et al. 2010, ’’Bailey et 
al. 2012, ’’Simon-Diaz et al. 2006, “’Stahl et al. 2008, ”Wade et al. 2006, ’^Grunhut et al. 2009, ’^Grunhut et al. 2012, ’'’Wade et al. 2012 


Table 2. GMRT observations of magnetic B and O type stars. 


Star 

Obs Date 

1.4 GHz 

Mean 

HJD 

Mean 

Phase 

Flux density 
F 20 (mJy) 

rms 

fiiy 

Obs Date 

0.6 GHz 

Mean 

HJD 

Mean 

Phase 

Flux density 
+50 (mJy) 

rms 

fiiy 

HD2I544I 

Ol-Nov-13 

2456598.262± 0.050 

0.93 

1.49*0.10 

46 

24-Oct-13 

2456590.168*0.044 

0.08 

0.98 *0.10 

56 

HD 37479 

31-Oct-13 

2456597.433* 0.170 

0.09 

2.00*0.11 

57 

24-Oct-13 

2456590.402* 0.044 

0.18 

1.25 *0.28 

51 

HD 37017 

31-Oct-13 

2456597.422* 0.131 

0.94 

1.59*0.12 

53 

24-Oct-13 

2456590.430* 0.040 

0.18 

0.59*0.32 

153 

HD 36485 

27-Oct-13 

2456593.378*0.050 

0.94 

1.30*0.10 

54 

26-Oct-13 

2456591.494* 0.076 

0.66 

0.59*0.17 

88 

HD 133880 

24-Jan-14 

2456681.616* 0.055 

0.84 

14.45*0.11 

49 

14-Jan-14 

2456671.590* 0.048 

0.42 

2.04*0.15 

66 

HD 37022 

31-Oct-13 

2456597.476 

0.4 

<4.95 

1650 






HD 57682 

Ol-Nov-13 

2456598.370 

0.13 

< 0.16 

52 






NGC 1624-2 

02-NOV-13 

2456598.487 

1.0 

< 0.21 

71 







The map resolutions for the 1390 MHz images were around 2”-3”, 
and around 5”-6” for the 610 MHz data. In Tabled we report the 
values of the final flux density. For non-detections, we quote 3(t 
upper limits. 


3 ANALYSIS AND RESULTS 

With our GMRT observations, we have detected 5 out of 8 stars. 
The detected stars, HD 215441, HD 37479, HD 37017, HD 36485 
and HD 133880, are B-type stars and have been detected in both the 
1390 and 610 MHz bands. The O-type stars, HD 37022, HD 57682 
and NGC 1624-2, were observed only in the 1390 MHz band, and 
none of them are detected. 

Of the B stars, 4 stars (HD 37479, HD 37017, HD 36485 and 
HD 133880) are sufficiently rapid rotators that rotation contributes 
significantly to the support of their magneto spheres (so-called ‘cen¬ 
trifugal magnetospheres’; [Petit et a02013ll . On the other hand, the 
O-type stars, as well as HD 215441, rotate slowly, and their mag¬ 
netospheres receive no sig nificant rotational support (so-called ‘dy¬ 
namical magnetospheres’; IPetit et aDl2013h . 

The observations are summarized in Table In this section, 
we give details of the detected B-type stars, and explore potential 
causes of non-detections in O-type stars. 


3.1 HD 215441 

HD 215441 (Babcock’s star) is a cool magnetic B8-9p star. It 
exhibits no optical, UV or X-ray evidence of a magnetosphere, 
but was det ected in the 5 GHz and 1.4 GHz bands (flux density 
~ 1 mjy) bv iLinskv et al.l ( Il992h . Despite its low temperature, this 


star is notewort hy due to its extraordinarily strong magnetic field 
(lBabcocklll96(]l) . 

Magnetic O- and B-type stars show variability of many obser¬ 
vational quantities. This variability is understood as rotation al mod¬ 
ulation in the context o f the oblique rotator model (ORM) dStibbj 
ll95d ; IWade et al.l201 ih . in which the dipolar surface magnetic field 
is tilted with respect to the rotational axis. Rotation of the star re¬ 
sults in strict periodic variability according to the stellar rotational 
period. 

To pha se our observations of HD 215441, we used the 
ephemeris of iNorth & Adelrn^ l ll995l) : 

ID = 2448733.714 -f (9.487574 ± 0.000030) • E. (1) 

Radio emission is clearly detected in both the 1390 and 610 
MHz bands at phases 0.93 and 0.08, respectively, with a flux den- 
sity somewhat larg er than that reported at higher frequencies by 
iLinskv et~^ ( Il992 . 1 .2-1.3 mJy). In both bands, our observations 
cover less than 1% of the rotational period; hence we are unable 
to evaluate any variability. The flux is somewhat (~ 15%) lower in 
the 610 MHz band. Whether this is due to the slope of the radio 
SED, or due to flux variability between phases 0.93 and 0.08, is not 
curr ently known. _ 

iLinskv et alJ i ll9921) reported a flux of 1.07 ± 0.09 mJy at 1.4 
GHz on July 2, 1987. According to Eq.[T] this corresponds to phase 
0.08, i.e. the same phase as our 610 MHz observation (0.98 ± 
0.1 mJy). Comparing the two measurements, this implies an essen¬ 
tially flat spectral index of (a = log(F’i 39 o/F 5 io)]/ log(1390/610) = 
+0.02 ± 0.03) at this phase. T his result appears i nconsistent with 
the a = -0.26 index derived bv iLeone et alj ) l2004l) based on higher 
frequency observations. This suggests that the SED flattens at low 
frequency. 





































Phase 


Figure 1. GMRT and VLA 1.4 GHz measurements (black and red dia¬ 
monds, respectively) and GMRT 610 MHz measurements (blue squares) of 
HD 37479 (cr Ori E). The VLA measurements are taken from the archive. 


3.2 HD 37479 


3.3 HD 37017 


HD 37017 is an early-type magnetic B2 Vp star that also exhibits 
Ho- emi ssion due to the p r esence of a dens e centrifugal magneto¬ 
sphere jLeone et al.ll201ol : IPetit et al.ll2Q13h . This target is in fact 
a double-lined spectroscopic binary (SB2) system with an 18.65d 
orbital period containing the magnetic B star in addition to a late B 
dwarf I Leoneet alj2010l). 

iLeone & Umanal 1 1993h detected variable radio emission from 
HD 37017. Like cr Ori E, the radio emission was found to vary with 
the rotational period, and m agnetic extrema coi ncided with radio 
maxima. It was detected by iLinskv et alj ( Il992h at 15, 5 and 1.4 
GHz bands using the VLA, and observed at several phases. The 
flux density was observed to vary by a factor of about 2, from 0.9- 
2.1 mJy at 1 5 GHz, 1-4-2 6 mjy at 5 GHz, and L5-2.4 mjy at 1.4 
GHz bands. iLeone et"^ l l2004l) studied HD 37017 in mm bands 
and examined the properties of the radio spectrum in the 1.4-87.7 
GHz range. They found the flux density to increase up to a fre¬ 
quency of 22.5 GHz and then to decrease in the mm range. This 
indicates a possible cut-off frequency. 

We obtained GMRT observations in the 1390 a nd 610 MHz 
bands . We phased the data using the ephemeris of iBolton et al.l 
l ll998l) : 


JD = 2446010.3750 -f (0.9011836 ± 0.000006) ■ E 


(3) 


HD 37479 (cr Ori E) is the prototypical B2 Vp star hosting a cen¬ 
trifugal magnetosphere. It exhibits strong an d variable Ho- emission 
that has been stud ied in some detail ( e.g. Townsend et '^ 120051 : 
lOksala et~^l2012h . It was detected bv iLinskv et al.l (Il992h at 15, 
5 and 1.4 GHz frequency bands using the VLA, and observed at 
several phases. The flux density was observed to vary by a factor 
of about 2, from 2.4-3.1 mJy at 15 GHz, 2.8-3.9 mJy at 5 GHz, and 
1.5-3.2 mJy at 1.4 GHz bands. They also detected variable circular 
polarization at 15 GHz. _ 

This star was also observed bv ILeone & IJmati3 lll993l) . who 
detected 5 GHz emission varying with the rotational period, with 
maxima coinciding with the magnetic extrema. 

We obtained multiple GMRT observ ations in both bands. T o 
phase the data, we use the ephemeris of iTownsend et al] l l2010l) . 
which takes into account the (slow) rotational braking of the star: 


JD = 2442778.8290 -f 1.1908229 • £ -H 1.44 x 10^^ • EL 


( 2 ) 


Our observations of HD 37479 were obtained at phases 0.95, 
0.08, 0.13 and 0.22 in the 1390 MHz band, and phases 0.15 and 
0.21 in the 610 MHz band. We have supplemented our observa¬ 
tions w ith the four 1.4 GHz measurements obtained bv iLinskv et al.l 
1 I 1992 I) . The phased flux variation is shown in Fig.[T] 

The combined 1390 MHz data suggest a roughly sinusoidal 
variation peaking at phase 0.5, with a peak flux of ~ 3.2 mJy. The 
610 MHz measurements are both weaker than the minimum 1390 
MHz measurement, but their phase coverage is too sparse to infer 
mea ningful information abou t the phase variation in this band. 

ILeone & Umanal d 19931) proposed that the 5 GHz flux from 
this star exhibited a double-wave phase variation. This is not sug¬ 
gested to be the case at 1.4 GHz band. Nevertheless, additional data 
will be required to draw more robust conclusions. 

Two of the 1390/610 MHz GMRT measurements were ob¬ 
tained at similar phases. From these we infer an instantaneous (0 = 
0.15-0.2) index of the SED at these frequencie s of or = 0.11 ±0.04 . 
This is consistent with a = 0.12 derived bv ILeone et al.l ( |2004|) 
based on higher frequency observations. 


Our 1390 MHz observations were obtained at phases 0.07, 
0.80 and 0.98, while our 610 MHz observations were obtained at 
phases 0.15 and 0.21. We suppl emented these measu rements with 
the 1420 MHz observations of iLinskv et"^ ( Il992l) . The phased 
data are illustrated in Fig.[^ 

ILeone & Umanal d 19931) described the 5 GHz flux variation as 
sinusoidal, with a maximum at phase 0.0 and a minimum at phase 
0.5. The variation shown in Fig. |2] although weak and coarsely 
sampled, appears opposite to this description. This is entirely a con¬ 
sequence of the different ephemeris used by those authors and by 
us. The phasing relative to the longitudinal field variation is the 
same. 

Again, our two 610 MHz observations exhibit lower flux than 
the lowest of the 7 measurements at 1390 MHz. One of the VLA 1.4 
GHz observations was obtained at a phase similar to the GMRT 610 
MHz data. From these we infer an instantaneous (0 = 0.15 - 0.2) 
index of the SED at these frequenci es of o- = 0.23 ±0 .1. This is 
consistent with a = 0.15 derived bv ILeone et al.l d2004h based on 
higher frequency observations. 

3.4 HD 36485 

HD 36485 {S Ori C) is a another hot B3 Vp star hosting a centrifu¬ 
gal magnetosphere. The Ha emission is variable on a timescale of 
a few hours. Notwithstanding its sharp spectral lines, the star is a 
rapid rotator with a period of about L5d, implying that it is viewed 
close to the rotational pole dBohlender et al.l 1 199 ll : ILeone et'H] 

l20l^ . 

6 Ori C was detected with the VLA at 5 GHz with a mean 
flux density of 0.95 mJy and showed a negative spectral index 
(a < -0.7; f y PC W ) between 5 and 15 GHz bands according 
to brake et"^ ) ll987l) . who concluded a non-thermal origin of the 
radio emission. Their 3 measurements at this wavelength did not 
show any significant variability, with a measurement uncertainty of 
±0.1 mJy. 

We use the ephemeris of iLeone et alj )l2010l) : 

JD = 2448298.86 ± (1.47775 ± 0.00004) • E. 


(4) 
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♦ GMRT1.4GHZ 

♦ VLA1.4GHZ 
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Figure 2. GMRT and VLA 1.4 MHz measurements (black and red dia¬ 
monds, respectively) and GMRT 610 MHz measurements (blue squares) of 
HD 37017. The VLA measurements are taken from the archive. 

In our GMRT observations, we detected this star in both the 
1390 and 610 MHz bands. In our 3 hour observation at 1390 MHz 
covering phases 0.91-0.96, the flux density varied between 0.42 ± 
0.24 to 0.6 ± 0.14 mjy. Although these measurements are formally 
consistent with a constant flux density, the flux did monotonically 
increase during the observation, suggesting a real change. In the 
610 MHz band, the flux density varied significantly, from 1.46 ± 
0.16 to 0.69 ± 0.13 mJy between phases 0.62 to 0.71. 

3.5 HD 133880 

HD 133880 (HR Lup) is a cooler, rapidly rotating {P = 0.88 d; 
ISailev et alj2012h Bp star hosting one of the strongest known stel¬ 
lar magnetic fields. However, unlike most Ap/Bp stars, HD 133880 
has a magnetic held topology that appears to be p r edominantly 
quadr upolar as opposed to dipolar ( iLandstreetll 199d : [Bailey et al.l 
I 2 OI 2 I) . HD 133880 was previously observed with the Australia 
Telescope Compact Array (ATCA) at 5 and 8.5 GHz frequency 
bands simultaneous ly on 1995 Febru ary 12, 14 and 16, for around 
10 hours each day. iLim et af] i ll9961) had demonstrated that radio 
flux and circular polarization of HD 133880 at both frequencies 
vary significantly and coherently according to the rotational pe¬ 
riod. They reported that the emission shows broad peaks near the 
ph ases of the lo n gitudi nal field extrema. The reanalysis of the data 
by iBailev et al.l ( I2012I) yielded the same result. They found flux 
variation to be complex, characterized by strong, broad maxima at 
phases 0.0 and 0.5 (i.e. the extrema of the longitudinal field), and 
sharper, somewhat weaker secondary extrema at quadrature phases 
(i.e. 0.25 and 0.75). _ 

ICeorge & Stevenj ( 12012l) observed this star with the GMRT in 
dual mode in the 610 and 240 MHz bands on December 5 and 7th, 
2009. They reported non-detections in all the observations. How¬ 
ever, we extracted the data from the GMRT archive and reanalysed 
it, and find detections on both days in the 610 MHz band. Further¬ 
more, we also find variability in each observation in this band. To 
re-confirm our detections of their data, we have analysed the data 
both in AIPS, as well as in Common Astronomy Software Appli¬ 
cations (CASA), separating out LL and RR Stokes, and we find 
consistent results. However, in the 240 MHz ba nd we obtained no 
detect ion in either observation, consistent with lGeorge & Stevens! 
( I 2 OI 2 I) . The 3 - cr upper limits in this band at the two epochs were 

1.5 and 3.0 mJy, respectively. We also located three data sets at 


1420 MHz in the VLA archive, obtained between 15-17 February 
1995. We reanalyse these data and tabulate them in Table [3 along 
with the flux density in various time intervals in both our new ob¬ 
servations, as well as the reanalysed observations. _ 

T o phase the data, we used the ephemeris of iBailev et al.l 

( I2OI2I) : 

ID = 2445472.0 -f (0.877476 ± 0.000009) • E. (5) 

The GMRT observations were made in the 1390 and 610 MHz 
bands. In Fig. 0 we plot the phased data in both bands. The flux 
variation of HD 133880 at both 1390 and 610 MHz is extraordi¬ 
nary. Unlike the other stars observed in this program (Table |4j that 
exhibit a flux variation of about a factor of 2, the 610 MHz flux of 
HD 133880 varies by a factor of 16 in our observations (and more 
than a factor of 12 at 1390 MHz). Even with the sparse phase cover¬ 
age of our observations, it is clear that the maxima of both the 610 
and 1390 MHz data occur at phases ~ 0.25 and 0.75, coincident 
with the phases of the minor maxima at 5 and 8 GHz described by 
IBailev etalJ ( l2012l) . 

Considering the important differences in the phase variations 
of the GMRT data and the ATCA data, we speculate that HD 
133880 may exhibit maser emission in a manner similar to CU Vir 
jTrigilio et al.ll20()(]|) . In CU Vir, while no coherent emission was 
found at 5 GHz, the 1.4 GHz emission was identified as ECME 
with a basal flux of 2-3 mJy and then very large increments in the 
flux density at around phases 0.35-0.45 and 0.75-0.85. While, we 
do not yet have complete phase coverage for HD 133880, the ob¬ 
served features are highly reminiscent of those of CU Vir (see Fig. 
0 . 

3.6 Free-free absorption in O star winds 

ISchnerr et alJ ( l2007l) carried out 5 and 1.4 GHz observations of 5 
(non-magnetic) O stars with the Westerbork Synthesis Radio Tele¬ 
scope (WSRT). They detected 3 stars, for the first time at 1.4 GHz 
band, with one star (f Per) inferred to show non-thermal radio emis¬ 
sion while two other (a Cam and A Cep) were inferred to show ther¬ 
mal radio emission. In all cases, the observed flux was lower than 
their predicted thermal flux, based on the mass loss rates inferred 
from the Her line profiles. They explained this discrepancy as due 
to a higher influence of wind clumping on the formation of the Her 
line. 

In our program we also observed 3 magnetic O-type stars at 
1390 MHz. None of these targets were detected (see Table0. For 
one of these stars (HD 37022) our upper limit is quite poor because 
it resides in a region with extended radio emission in Orion. 

The high mass loss rates of O stars can make the wind op¬ 
tically thick due to free-free absorption out to a very large “radio 
photosphere”. In such a case radio observations are unlikely to di¬ 
rectly see the magnetospheric emission. On the other hand, they 
may probe the global mass loss, which is predicted to be quite 
heavily reduced for O stars with dynamical magnetospheres, due 
to que nching of the mass loss by the field and the result ant plasma 
in-fall dud-Doula & OwocOl2002l : [ud-Doula et al.ll2008l) . 

We can estimate the radius of the free-free radio photosphere, 
under the assumptio n of a spheric ally-symmetric, non-magnetic 
wind, using Eq. (4) of lTorreslO^ : 

Tff = 5 X (6) 

where Tg is the free-free optical depth, M is the mass-loss rate in 
units of 10“* 3fg/yr, Voo is the wind terminal velocity in units of 
10* cm/s, / is the clump volume filling factor, v is the frequency of 
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Table 3. Variability of HD 133880 . 


Freq 

GHz 

Telescope 

Date 

Mean 

HJD 

Phase 

Flux density 
mJy 

rms 

fiiy 

1.388 

GMRT 

24-Jan-14 

2456681.572 ±0.012 

0.791 

25.90 ±0.14 

108 



24-Jan-14 

2456681.604 ±0.012 

0.827 

26.14 ±0.12 

123 



24-Jan-14 

2456681.636 ±0.012 

0.863 

3.41 ±0.12 

117 



24-Jan-14 

2456681.663 ±0.008 

0.895 

2.43 ±0.14 

177 

1.425 

VLA 

I5-Feb-95 

2449764.040 ± 0.007 

0.348 

2.07 ± 0.40 

396 



15-Feb-95 

2449764.098 ± 0.006 

0.414 

2.87 ± 0.33 

331 



16-Feb-95 

2449765.013 ± 0.006 

0.456 

3.39 ± 0.24 

239 



I7-Feb-95 

2449766.010 ± 0.006 

0.593 

5.86 ± 0.27 

267 

0.608 

GMRT 

14-Jan-14 

2456671.552 ±0.010 

0.371 

1.91 ±0.23 

76 



14-Jan-14 

2456671.579 ±0.010 

0.402 

1.97 ±0.24 

68 



14-Jan-14 

2456671.607 ±0.010 

0.434 

2.21 ± 0.24 

65 



14-Jan-14 

2456671.631 ±0.007 

0.461 

1.91 ±0.14 

75 

0.606 

GMRT 

05-Dec-09 

2455170.712 ±0.026 

0.966 

2.37 ±0.31 

311 



05-Dec-09 

2455170.776 ±0.027 

0.038 

2.16 ±0.33 

328 

0.607 

GMRT 

07-Dec-09 

2455172.663 ±0.027 

0.189 

1.26 ±0.32 

218 



07-Dec-09 

2455172.726 ±0.026 

0.261 

20.08 ± 0.40 

287 



07-Dec-09 

2455172.787 ±0.024 

0.330 

3.52 ± 0.37 

347 


Table 4. Variability of other stars in our sample. 


Star 

Tel. 

Freq 

GHz 

Date 

Mean 

HJD 

Phase 

Flux density 
mJy 

rms 

j'Jy 

Ref. 

HD215441 

GMRT 

1.388 

01 Nov 13 

2456598.223 ± 0.010 

0.927 

1.43 ±0.18 

76 

This paper 




01 Nov 13 

2456598.255 ± 0.010 

0.931 

1.15±0.15 

70 





01 Nov 13 

2456598.283 ± 0.010 

0.934 

1.45 ±0.18 

74 





01 Nov 13 

2456598.306 ± 0.005 

0.936 

1.20 ±0.19 

92 



VLA 

1.4 

02 Jul 87 

2446979.250 

0.078 

1.07 ±0.09 



HD215441 

GMRT 

0.608 

24 Oct 13 

2456590.136 ±0.012 

0.075 

1.19±0.13 

71 

This paper 




24 Oct 13 

2456590.168 ±0.012 

0.078 

1.18±0.13 

69 





24 Oct 13 

2456590.201 ± 0.012 

0.082 

0.70 ± 0.22 

108 


HD37479 

GMRT 

1.388 

31 Oct 13 

2456597.273 ± 0.010 

0.951 

1.45 ±0.17 

87 

This paper 




31 Oct 13 

2456597.430 ± 0.010 

0.082 

1.62 ±0.15 

84 





01 Nov 13 

2456597.514 ±0.010 

0.153 

2.01 ±0.17 

86 





01 Nov 13 

2456597.595 ± 0.008 

0.221 

1.91 ±0.19 

106 



VLA 

1.4 

11 Mar 85 

2446136.467 

0.585 

3.2 ±0.4 


Linskv et al. G992) 




12 Mar 85 

2446136.685 

0.768 

2.7 ±0.3 






16 Mar 85 

2446141.480 

0.795 

2.1 ±0.3 






17 Mar 85 

2446141.699 

0.979 

1.5 ±0.3 



HD37479 

GMRT 

0.608 

24 Oct 13 

2456590.370 ± 0.012 

0.154 

1.25 ±0.24 

134 

This paper 




24 Oct 13 

2456590.434 ± 0.012 

0.208 

0.80 ±0.33 

158 


HD37017 

GMRT 

1.388 

31 Oct 13 

2456597.301 ± 0.010 

0.802 

1.08 ±0.15 

74 

This paper 




31 Oct 13 

2456597.458 ± 0.010 

0.975 

1.04 ±0.13 

66 





01 Nov 13 

2456597.542 ± 0.010 

0.069 

1.54 ±0.16 

78 



VLA 

1.4 

11 Mar 85 

2456136.488 

0.459 

2.4 ± 0.2 


Linskv et al. (1992) 




12 Mar 85 

2456136.623 

0.609 

1.8 ±0.2 






16 Mar 85 

2456141.419 

0.931 

1.5 ±0.3 






17 Mar 85 

2456141.637 

0.173 

1.5 ±0.3 



HD37017 

GMRT 

0.608 

24 Oct 13 

2456590.402 ± 0.012 

0.146 

0.50 ±0.37 

74 

This paper 




24 Oct 13 

2456590.462 ± 0.008 

0.212 

0.42 ± 0.24 

66 


HD36485 

GMRT 

1.388 

27 Oct 13 

2456593.338 ± 0.010 

0.910 

0.43 ±0.17 

90 





27 Oct 13 

2456593.365 ± 0.010 

0.929 

0.44 ± 0.33 

85 





27 Oct 13 

2456593.393 ± 0.010 

0.947 

0.60 ±0.14 

92 





27 Oct 13 

2456593.419 ± 0.009 

0.965 

0.52 ±0.17 

78 



VLA 

1.4 

07 Jul 86 

2446619.810 

0.221 

< 1.2 

400 

Linskv et al. G992) 

HD36485 

GMRT 

0.608 

25 Oct 13 

2456591.429 ±0.012 

0.618 

1.46 ±0.16 

108 

This paper 




25 Oct 13 

2456591.461 ±0.012 

0.640 

1.26 ±0.18 

218 





26 Oct 13 

2456591.526 ±0.012 

0.684 

0.96 ±0.19 

87 





26 Oct 13 

2456591.558 ±0.012 

0.706 

0.69 ±0.13 

97 
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Figure 3. Upper frame - Phase versus flux density plot of HD 133880 at 1.4 
GHz (red color) and 610 MHz (blue color). Here red squares denote GMRT 
1390 MHz measurements, red diamonds denote VLA 1.4 GHz measure¬ 
ments, blue squares denote GMRT 610 MHz observations made by us in 
2014 January, whereas th e blue diamonds are the 2009 December data of 
iGeorge & Stevenj j2012h reanalyzed by us. Note that unlike in Figs. 1 and 
2, the vertical scale in this figure is logarithmic. Lower frame - Long itudi- 
nal magnetic field variation of HD 133880, from iBailev et alj feoi j) . The 
vertical lines in both plots indicate phases 0.25 and 0.75. 


observation in GHz, T is the wind temperature in units of 10^ K, 
which we here for simplicity approximate with the effective tem¬ 
perature of the star, and Rg is the distance from the star, in units of 
3 X 10*^ cm. 

We adopt a volume filling factor / = 1, i.e. a smooth, un- 
clumped radio e mitting region, since it has been demonstrated 
JPuls et alj I2OO6I) that f = I yields the best agreement with 
observations of thermal r a dio em ission from hot stars when the 
IVink, de Koter, & Lamerd ( I2OOOI) scaling law for mass-loss rates 
is assumed. Assuming Tff = 1, we have computed the position of 
the radio photosphere (£), in units of Rf) at 1390 MHz, using the 
wind parameters (com puted using the V ink scaling) and physical 
parameters reported by [Petit et al J ( 12013), except for HD 133880, 
for which we use values reported by [Bailey et al.l ( l2012h . The re¬ 
sults are summarized in Table[T]and illustrated in Fig. [4] 

If Rfs < Ra we expect the magnetospheric emission to escape 
the wind and to be detectable (in principle). However, if Rff > Ra, 
the magnetospheric emission will be hidden within the (thermal) 
radio photosphere of the wind. 

In our sample, B stars show ratios of Rb/Ra which are clus¬ 
tered around values 0.1-1.5, whereas all the O-type stars have this 
ratio clustered between 10-100. This is illustrated in Fig. |4] This 
convincingly demonstrates that the free-free absorption of the mag¬ 
netospheric radio emission of the B stars likely is negligible, and 


that the principal cause of the null results for O stars is free-free 
scattering in the free streaming wind. 

Current models of the magnetically confined winds of slowly 
rotating O stars imply that they host “dynamical magnetospheres”. 
In this scenario, wind streamlines inside Ra in the magnetic equa¬ 
torial region are constrained to follow closed field loops. After the 
wind stream reach the respective loop top and shock, the cooling 
plasma returns to the stella r surface on a relatively sho rt (dynami¬ 
cal) timescale. Simulations dud-Doula & Owockil2002l) predict that 
as a consequence of this confinement and fail-back, the global mass 
loss of O stars with dynamical magnetospheres should be signif¬ 
icantly reduced, by up to a factor of 10-20 compared to similar 
non-magnetic stars. 

In the context of free-free absorption, this reduction in mass 
loss has important implications for the location of the radio photo¬ 
sphere. T he expected mass-loss rate can be estimated using equa¬ 
tion 23 of lud-Doula et alj (l2008h . leading to a reduction of M by a 
factor of more than 3 for HD 37022, 5 for HD 57682, and 15 for 
NGC 1624-2. If we reduce the mass loss rate in Eq. (6) by these 
factors, for HD 37022 and HD 57682 the ratio Rb/Ra is reduced to 
roughly 40, i.e. still well outside the Alfven radius. However, the ra¬ 
tio for NGC 1624-2 is reduced to ~ 2. At shorter wavelengths (3-6 
cm), this value becomes < 1. Hence multiwavelength radio obser¬ 
vations spanning the radio spectrum have the potential to provide 
novel constraints on predictions of mass-loss quenching of strongly 
magnetic O stars like NGC 1624-2. 


4 SUMMARY AND DISCUSSION 

We have studied 8 O- and B-type stars with the GMRT in the 1390 
and 610 MHz bands. Our 610 MHz detections of the 5 B type stars 
are the first detections of this class of objects at such low frequen¬ 
cies. This indicates that free-free absorption of their winds is not 
able to shadow their low frequency emission. We also reanalysed 
archival o bservations of HD 133880 obtained in December 2009, 
for which iGeorge & StevensI ( 1201 2h claimed non-detections. Sur¬ 
prisingly, we do detect the radio emission in both their observations 
allowing us to diagnose variability. 

For HD 37017 and HD 37479, we confirm weak variab ility of 
the 1390 MHz flux density, as reported bv lUeone & Umanj ( Il993h 
in 5 GHz band. The variability appears to be in phase with the vari¬ 
ation of the longitudinal magnetic field. We do not find evidence of 
a double-wave 1390 M Hz variation of HD 37479 (as reported by 
lUeone & Umanairi993L at 5 GHz), although better phase coverage 
is necessary to robustly confirm this conclusion. 

For the B star HD 36485, our observations suggest a relatively 
rapid change in the flux density in both bands in narrow phase 
ranges. This result needs to be confirmed by further observations. 

The cooler B star HD 133880 shows remarkable variability in 
both the 1390 and 610 MHz bands. Even with our limited phase 
coverage, the initial results suggest that the fluxes in both bands 
peak at pha ses coinciding with the minor peaks obtained from 
ATCA data dBailev et alj]2012h . We propose that this behaviour 
may be a consequence of ECM E, first observed from the B8p star 
CU Vir bv lTrigilio et al.l ( l2000h . The phenomenon, only observable 
at low radio frequencies, strongly motivates immediate monitoring 
of this target. 

There were no detections obtained for any of the thre e mag¬ 
netic O stars in the 1390 MHz band. Using the scaling law of lTorrej 
( l201lh . we are able to conclude that this null result is in gen¬ 
eral agreement with the expected free-free optical depth of dense 
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Figure 4. Ratio of Rg to Ra versus effective temperature for the observed 
sample (excluding HD 215441, for which Ra is unavailable). Here red cir¬ 
cles and blue circles are for B- and O- type stars. The inverted triangle 
indicate upper limits as listed in Table[T] 



winds of 0-type stars. We also conclude that radio observations at 
higher frequencies will be capable of testing predictions of mass 
loss quenching of strong ly magnetic O stars like NGC 1624-2. 

iLeone et aP ||2004 ) commented that the radio spectra of HD 
142301 and HD 215441 (the two stars in their sample with the 
strongest magnetic fields) exhibited a sharp drop in flux density 
between 15 and 22.5 GHz. They concluded that this suggests there 
is a cut-off frequency in this range. On the other hand, the mag¬ 
netic chemically peculiar stars with relatively weak magnetic fields 
showed flux densifies fhat cerfainly decrease from the cm to the mm 
range. In our sample, we were able to measure spectral indices be- 
tween 1390 and 61 0 MHz for 3 stars in common with the sample of 
iLeone et~^ I ( |2004 . For two stars (HD 37017 and HD 37479), the 
indices were similar to those obtained at higher frequencies, sug¬ 
gesting a smooth decrement of the flux density. On the other hand, 
for HD 215441 (the nondegenerate star with the strongest known 
magnetic field), we o btained a spectral in dex that differed strongly 
from that derived by ILeone et"^ (l2004h . This implies an equally 
significant change in the spectral index at low frequencies, and de¬ 
serves further investigation. 

The results presented in this paper represent a first feasibility 
study of the detection of low frequency emission from magnetic 
O- and B-type stars. Having demonstrated the general properties 
of a range of objects, we are now carrying out a larger survey to 
observe all known magnetic OB stars in radio bands from low to 
high frequencies. We are also performing monitoring of individual 
magnetic stars, in order to better understand the rotational variation 
of their fluxes and spectral energy distributions. 
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